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ABSTRACT 

Galaxy groups are not scaled down versions of massive galaxy clusters - the hot 
gas in groups (known as the intragroup medium, IGrM for short) is, on average, less 
dense than the intracluster medium, implying that one or more non-gravitational pro- 
cesses (e.g., radiative cooling, star formation, and/or feedback) has had a relatively 
larger effect on groups. In the present study, we compare a number of cosmological 
hydrodynamic simulations that form part of the Overwhelmingly Large Simulations 
project to isolate and quantify the effects of cooling and feedback from supernovae 
(SNe) and active galactic nuclei (AGN) on the gas. This is achieved by comparing La- 
grangian thermal histories of the gas in the different runs, which were all started from 
identical initial conditions. While radiative cooling, star formation, and SN feedback 
are all necessary ingredients, only runs that also include AGN feedback are able to 
successfully reproduce the optical and X-ray properties of groups and low-mass clus- 
ters. We isolate how, when, and exactly what gas is heated by AGN. Interestingly, we 
find that the gas that constitutes the present-day IGrM is that which was not strongly 
heated by AGN. Instead, the low median density/high median entropy of the gas in 
present-day groups is achieved by the ejection of lower entropy gas from low-mass 
progenitor galaxies at high redshift (primarily 2 < z < 4). This corresponds to the 
epoch when supermassive black holes accreted most of their mass, typically at a rate 
that is close to the Eddington limit (i.e., when the black holes are in a 'quasar mode'). 

Key words: galaxies: formation — galaxies: clusters: general — galaxies: groups: 
general cosmology: theory X-rays: galaxies: clusters intergalactic medium 



1 INTRODUCTION 

Understanding the formation and evolution of galaxy groups 
is crucial to solving the general problem of galaxy forma- 
tion, as groups contain most of the galaxies in the universe 
at the present day (e.g., Mulchaey 2000), and a large frac- 
tion of the universal stellar mass is formed in groups (e.g., 
Crain et al. 2009). They can also potentially yield important 
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insights into the formation and evolution of more massive 
and rare galaxy clusters, since in the current cosmological 
paradigm clusters acquire most of their mass through the 
accretion/mergers of groups. Groups of galaxies are special 
systems in an observational sense as well, as it is possible 
to probe both the stellar and hot gas components of these 
systems directly, whereas for normal galaxies detecting the 
diffuse hot gas is presently very challenging. As highlighted 
in many recent theoretical papers (e.g., Bower et al. 2008, 
hereafter B08; McCarthy et al. 2010; Fabjan et al. 2010; 
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Puchwein et al. 2010), the properties of the hot gas and of 
the stellar component are intimately linked and any success- 
ful model must be able to reproduce both simultaneously. 

X-ray observations of nearby, relatively bright galaxy 
groups have revealed that the hot gas (the intragroup 
medium, IGrM for short) has some puzzling properties. In 
particular, it is on average less dense than the hot gas found 
in massive clusters (known as the intracluster medium, ICM 
for short). This is not what one expects if gravitational pro- 
cesses alone determine the thermodynamic properties of the 
gas (e.g., Voit, Kay, Bryan 2005, hereafter VKB05). The ob- 
served decrease of the gas density with halo mass therefore 
indicates that one or more non-gravitational processes (e.g., 
radiative cooling, star formation, and/or feedback) has had 
a larger effect on groups than on massive clusters. The fact 
that a non- negligible fraction of the total mass of baryons 
in groups is locked up in stars (e.g., Balogh et al. 2001; Lin 
& Mohr 2004; Gonzalez et al. 2007; Giodini et al. 2009) is 
already testament that non-gravitational processes are im- 
portant on the scale of groups. 

The entrop30 of the IGrM is perhaps the best tool for 
studying the impact of non-gravitational physics, since it 
maintains a record of the thermal history of the gas (see 
Voit et al. 2003 and Voit 2005 for further discussion). Using 
ROSAT observations, Ponman et al. (1999) were the first 
to show convincingly that galaxy groups have 'excess' en- 
tropy relative to massive clusters; i.e., when one scales the 
radius and entropy by some suitable characteristic radius 
and entropy that depend only on system mass (e.g., the 
virial radius and 'virial entropy'; see Section 3.1), groups 
typically have higher scaled entropies at all scaled radii 
than clusters. And both groups and clusters have excess en- 
tropy relative to theoretical models that only account for 
gravitational heating (e.g., Balogh et al. 1999; Babul et al. 
2002; Voit et al. 2002). More recent observations, partic- 
ularly with ASCA, Chandra, and XMM-Newton, have im- 
proved the quality of the measurements and extended both 
the range of the masses and radii over which the entropy 
has been measured (e.g., Finoguenov et al. 2002; Ponman et 
al. 2003; Sun et al. 2009; Johnson et al. 2009; Pratt et al. 
2010; Giodini et al. 2010). These newer results have generally 
confirmed those from the ROSAT satellite: groups contain 
excess entropy relative to clusters as well as to models that 
only take into account gravitational heating. 

On the theory side, many suggestions have been put for- 
ward for the physical origin of the excess entropy. The most 
commonly mentioned are 'preheating' (first introduced by 
Kaiser 1991 and Evrard & Henry 1991), radiative cooling 
(e.g., Bryan 2000), and supernova (SN) and active galactic 
nuclei (AGN) feedback within groups. Preheating is concep- 
tually simple: at early times the entropy of the proto-lGrM is 
raised by feedback processes (such as outflows driven by SN 
or AGN). The increased entropy prevents the gas from be- 
ing easily accreted or compressed, resulting in a lower mean 



1 As is common in X-ray astronomy, we define the 'entropy', S, 
as fcsT/ne^ 3 and use units of keV cm 2 , where ks is Boltzmann's 
constant, T is the temperature of the gas and n e is the electron 
number density. The observational entropy S is related to the 
thermodynamic specific entropy s via s oc In S and, like the ther- 
modynamic entropy, S will be conserved in any adiabatic process. 



density (and gas mass fraction) for groups. Note that for a 
fixed level of entropy injection the effect is larger for groups 
than for clusters, since groups have lower entropies than 
clusters. Somewhat paradoxically, radiative cooling has also 
been shown to raise the entropy of the gas in groups and 
clusters. It does so not by heating the gas, but by selec- 
tively removing the lowest entropy gas from the hot IGrM 
phase, turning it into cold gas and stars, and thereby rais- 
ing the mean entropy of the hot gas that is unable to cool 
(e.g., Bryan 2000). Feedback from winds driven by SNe and 
jets/winds launched by supermassive black holes (hereafter 
BHs) originating from orbiting group galaxies may also sig- 
nificantly raise the entropy of the IGrM by direct heating. 
At present, no clear consensus has emerged as to which of 
these processes are most important for raising the entropy of 
the gas over and above that generated through gravitational 
shock heating alone. 

In the present study, we seek to develop a better under- 
standing of the physical processes that set the global thermo- 
dynamic properties of the IGrM. We will do this by exploit- 
ing a suite of high resolution cosmological hydrodynamical 
simulations which form part of the Overwhelmingly Large 
Simulations project (OWLS; Schaye et al. 2010). The pur- 
pose of the OWLS project is to explore the effects of varying 
so-called 'subgrid' physics in the simulations, such as radia- 
tive cooling, star formation, and SN and BH feedback. All 
simulations are started from identical initial conditions, with 
each simulation varying a parameter within one of the sub- 
grid physics modules or switching on/off one of the subgrid 
processes entirely. By comparing the different simulations, 
it is possible to isolate which processes are most important 
for raising the entropy of the IGrM. In particular, we will 
compare not only the present-day entropy distributions of 
the groups in the different simulations, but also, for the first 
time, the Lagrangian histories of the gas in the different 
runs. This allows us to follow the heating and cooling of a 
parcel of gas over cosmic time. We will show that a partic- 
ularly powerful discriminator is to identify and follow the 
thermal properties of the same sets of particles in the (iden- 
tical) initial conditions of the various simulations. 

In McCarthy et al. (2010; hereafter Paper I) we showed 
that the OWLS AGN model yields an excellent match to 
a wide range of optical and X-ray observations of nearby 
galaxy groups, including the entropy, temperature, and 
metallicity distribution and their dependence on halo mass, 
as well as the stellar mass of the group and the central 
brightest galaxy (CBG) and the stellar age of the CBG. 
Based on comparisons between this model and a number 
of other models, we argue below that, in fact, none of the 
three physical scenarios outlined above provides an accurate 
description for how the entropy of the IGrM was raised. In- 
stead, we show that gas that constitutes the present-day 
IGrM was that which escaped significant AGN feedback and 
the excess entropy is due primarily to the ejection of the 
lowest entropy gas at high redshift (i.e., the gas that was 
significantly heated is no longer within the virialized region). 

The present paper is organised as follows. In Section 2 
we provide a brief description of the simulations. In Section 
3 we present our main results on the thermal history of the 
IGrM in the simulations. In Section 4 we explore further 
how the AGN feedback raises the entropy of the IGrM. Fi- 
nally, in Section 5 we summarize and discuss our findings. 
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All quantities presented in this study (e.g., entropy, tem- 
perature, density) are in physical, rather than co-moving, 
units. 



2 SIMULATIONS 

An in-depth presentation of the OWLS project can be found 
in Schaye et al. (2010). We therefore present only a brief 
description of the simulations here. 

All of the OWLS runs we analyse are initialised from 
identical initial conditions, in particular from a ACDM cos- 
mological density field (with a power spectrum computed 
using CMBFAST) in a periodic box of WOh' 1 Mpc on a 
side at z = 127. In generating this density field, the vari- 
ous relevant cosmological parameters (fit, f2 m , S7a, h, erg, 
n s ) were set equal to the maximum-likelihood values found 
from the analysis of the 3- year WMAP cosmic microwave 
background data (0.0418, 0.238, 0.762, 0.73, 0.74, 0.95, re- 
spectively; Spergel et al. 2007). The simulations are evolved 
to 2 = using the TreePM-SPH code GADGET-3 (last de- 
scribed in Springel 2005) . We extract from the final snapshot 
(z — 0) all galaxy groups for which log lo (M2oo/M0) ^ 13.00 
(where M200 is the total mass within a radius that encloses 
a mean density of 200 times the present-day critical den- 
sity of the universe). The simulations use 2 x 512 3 parti- 
cles, yielding a particle mass for gas and dark matter of 
m gas « 8.65 x 10 7 /i _1 M Q and m dm ~ 4.06 x 10 s /T 1 M Q , 
respectively. Thus, a ~ 10 14 ft -1 Mq galaxy group is resolved 
with ~ 10 5 gas and dark matter particles. As shown in the 
Appendix of Paper I, this resolution is sufficient to robustly 
predict the properties of galaxy groups in the simulations. 

In Paper I (see also Schaye et al. 2010), we described in 
some detail the implementation of 'subgrid' physics (e.g., ra- 
diative cooling, star formation, and SN and AGN feedback) 
in the GADGET-3 code and how we select and analyse galaxy 
groups. We refer the reader to that study for details. As dis- 
cussed in Section 1, the thermal properties of the IGrM differ 
from those predicted by simulations that include only grav- 
ity and basic non-radiative hydrodynamics. As a result, it 
is now widely believed that one or more non-gravitational 
processes are responsible for modifying the properties of the 
hot gas. The most likely culprits are cooling, star formation, 
and feedback, which are implemented in the OWLS simula- 
tions using subgrid prescriptions. The power of the OWLS 
project is not only that it systematically varies the relevant 
physical parameters of a given subgrid physics module, but 
that it also systematically varies which subgrid processes are 
switched on/off. This is the best and perhaps only way to 
isolate which processes are most important for establishing 
the thermodynamics properties of the IGrM. 

We note that the physics that raises the entropy of the 
ICM in massive clusters is quite likely to be closely related 
to that which operates on the scales of groups. However, the 
volume of the OWLS runs [(100ft -1 Mpc) 3 ] is not sufficiently 
large to contain massive galaxy clusters. We have plans to 
carry out both 'zoomed' initial condition and larger periodic 
box simulations in the future. 

Below, we provide brief descriptions of the runs that 
we analyse in the present study. The runs we have selected 
represent a range of interesting limiting cases and not the 
entire OWLS suite of simulations. In Table 1 we list all the 



simulations, the corresponding names in the OWLS project 
and the differences in the included subgrid physics. 

2.1 Individual runs 

• NoCool: This model omits the effects of radiative cool- 
ing, star formation, and feedback from SNe and BHs. It 
does, however, include net photoheating from a Haardt & 
Madau (2001) UV/X-ray background (i.e., whenever A nc t = 
Ahcat — A coo i > 0, otherwise A not = 0). However, this photo- 
heating has negligible consequences for the gas that consti- 
tutes the IGrM at the present-day for redshifts of z < 2 (see 
the Appendix). Therefore, this model is basically equivalent 
to a pure 'non-radiative' run below zr;2. We include net 
photoheating purely for consistency, since all other OWLS 
runs include UV photoheating. This model will serve as our 
baseline for quantifying how radiative cooling, star forma- 
tion, and feedback from SNe and AGN raise the entropy of 
the IGrM over and above that of gravitational shock heating. 

• PrimCooLSF: This model is identical to the ZC00LSF 
model (described below), except that cooling rates are cal- 
culated assuming primordial abundances. 

• ZC00LSF: Radiative cooling rates for the gas are com- 
puted on an element-by-element basis by interpolating 
within pre-computed tables (generated with the CLOUDY 
software package, last described in Ferland et al. 1998) that 
contain cooling rates as a function of density, tempera- 
ture, and redshift calculated in the presence of the cos- 
mic microwave background and the Haardt & Madau (2001) 
UV/X-ray background (see Wiersma et al. 2009a). Star for- 
mation is tracked in the simulations following the pressure- 
based prescription of Schaye &: Dalla Vecchia (2008), which 
reproduces the observed Kennicutt-Schmidt law by con- 
struction. The timed release of individual elements by both 
massive (Type II SNe and stellar winds) and intermedi- 
ate mass stars (Type la SNe and asymptotic giant branch 
stars) is included following the prescription of Wiersma et al. 
(2009b). A set of 11 individual elements are followed in these 
simulations (H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe), which 
represent all the important species for computing radiative 
cooling rates (Wiersma et al. 2009a). As is the case for other 
OWLS runs, UV/X-ray photoheating is included. However, 
as already noted, this is unimportant below z « 2 and this 
model therefore effectively represents a 'no feedback' model 
for galaxy groups. (Note that even though there is no form 
of effective feedback in this run, metals are still removed 
from galaxies and dispersed into the IGrM by ram pressure 
stripping.) 

• PrimCooLSFSN. This model is identical to the Prim- 
CooLSF model except that feedback from SN-driven winds 
is also included. In particular, the kinetic wind model of 
Dalla Vecchia & Schaye (2008) is adopted with the wind 
velocity, v w , set to 600 km/s and the mass- loading param- 
eter (i.e., the ratio of the mass of gas given a velocity kick 
to that turned into newly formed star particles), r\, set to 
2. This choice of parameters corresponds to using approxi- 
mately 40% of the total energy available SNe, for a Chabrier 
(2003) stellar initial mass function (IMF) (which is adopted 
in all of the simulations presented in this paper, with the 
obvious exception of NoCool). 

• ZCooLSF_SN. This model is identical to Prim- 
CoolJSFJSN except cooling rates are calculated on an 
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Table 1. A list of all simulations, the corresponding names in the OWLS project (Schaye et al. 2010) and the differences in the included 
subgrid physics. 



Simulation 


OWLS name 


Cooling 


Star formation 


Supernova 


AGN 










Feedback 


Feedback 


No Cool 




None 


No 


No 


No 


PrimCooLSF 


NOSN-NOZCOOL 


Primordial 


Yes 


No 


No 


ZCooLSF 


NOSN 


Metal 


Yes 


No 


No 


PrimCool.SF.SN 


NOZCOOL 


Primordial 


Yes 


Yes 


No 


ZCool.SF.SN 


REF 


Metal 


Yes 


Yes 


No 


ZCooLSF.SN.AGN 


AGN 


Metal 


Yes 


Yes 


Yes 



element-by-element basis, as described above for the 
ZCooLSF model. ZCool.SF.SN was referred to as REF in 
Paper I (see Table 1). 

• ZCool.SF.SN.AGN: This model is identical to the 
ZCool.SF.SN model except that it also includes the pre- 
scription for BH growth and AGN feedback of Booth & 
Schaye (2009), which is a modified version of the algorithm 
developed by Springel et al. (2005). A brief description is 
as follows. BHs grow in mass through mergers with other 
BHs and through the accretion of neighbouring gas. Gas is 
assumed to accrete at the local Bondi-Hoyle-Lyttleton rate. 
However, if the local gas density exceeds n H = 0.1 cm ~ 3 
then this rate is multiplied by the factor [riH/(0.1 cm~ 2 )] 2 
to compensate for the fact that our simulations lack both 
the resolution and the physics to model a cold, interstellar 
gas phase (see Booth & Schaye 2009 for a discussion). A 
certain fraction of the rest mass energy of the accreted gas, 
e, is used to heat a number of randomly selected neighbour- 
ing gas particles, Wheat, by raising their temperatures by an 
amount ATheat. (Note that the BHs store feedback energy 
until it is sufficient to heat Wheat particles by ATheat-) This 
efficiency is a free parameter of the model. Booth & Schaye 
(2009) find that a value of e = 0.015 yields a good match 
to the observed z — relations between BH mass and stel- 
lar mass and velocity dispersion and the z — cosmic BH 
density. We therefore fix e to this value. As discussed below 
(see Section 4.4), the predicted properties of galaxy groups 
are insensitive to the adopted value of riheat (which, for the 
present study, we set to unity) but are sensitive to ATheat if 
this quantity is set to a value that is comparable to (or lower 
than) the virial temperature of the high redshift progenitors 
of galaxy groups. For specificity, we set ATheat = 10 K and 

Wheat = 1 • 

To give a visual example of the importance of includ- 
ing non-gravitational processes in simulations of the IGrM 
we present Fig. 1, which shows projected mass- weighted 
entropy, surface gas mass density, and projected mass- 
weighted temperature of the gas in and around the same 
randomly-selected galaxy group in NoCool, ZCool.SF.SN, 
and ZCool.SF.SN.AGN. The maps arc 10 tr 1 Mpc on a 
side. Particles were interpolated to a 3D uniform grid (with 
1000 3 pixels) using the SPH smoothing kernel and projected 
along the z-axis of the simulation box. In NoCool the hot, 
high entropy gas is confined to the area in the immediate 
vicinity of the group and to a large filament that runs from 
the bottom left to the top right of the maps. The intro- 
duction of supernova feedback (and cooling and star forma- 
tion; ZCool.SF.SN) results in a relatively small amount of 



moderate-entropy gas being expelled into voids. The inclu- 
sion of AGN feedback has a much more dramatic effect, with 
large amounts of high entropy/high temperature gas being 
ejected out into the surrounding intergalactic medium. We 
discuss in detail the importance of gas expulsion by AGN in 
Sections 3 and 4. 



3 THE THERMAL HISTORY OF THE IGRM 

In this section we first present the radial entropy distribu- 
tions of the various OWLS runs, demonstrating that the in- 
clusion of cooling and feedback processes indeed raises the 
entropy of the hot gas in groups. We then examine the en- 
tropy histories of these runs in detail and compare them with 
one another to deduce which processes are most important 
for raising the entropy of the gas. 

3.1 Radial distribution of the entropy 

We begin by examining the 3D radial entropy distribu- 
tions of the hot gas (T > 10 s K) in the OWLS runs de- 
scribed in Section 2.1. In Fig. 2 we plot the median 3D 
mass-weighted entropy profiles for all groups with masses 
between 13.25 ^ log lo (Msoo/M0) ^ 14.25 [approximately 
70 groups with a median mass of log 10 (M5oo/M Q ) « 13.5], 
where M500 is the total mass within the radius that encloses 
a mean density of 500 times the present-day critical density 
of the universe. The radial coordinate has been normalised 
by r 5 oo, which is the radius that encloses M500 (typically 
rsoo ~ 0.65r2oo)- The entropy has been normalised by S500, 
the 'virial entropy', which we define, following VKB05, as: 

or n\ — fcsTioo , . 

= °> = i^cF 3 " (1) 

GM 500 ij,mH 
2r 5 oo[500 fbPcrit(z = 0)/(u e m H )] 2 / 3 



where ft, = f2j,/f2 m , and jj, e is the mean molecular weight per 
free electron. Note that S500 is not the entropy at rsoo, i.e., 
S500 7^ S(r 500), even in the case of the self-similar model. 
This is because S500 is defined in terms of the mean gas 
density within rsoo (as opposed to the gas density at rsoo) 
for a halo with the universal baryon fraction and the virial 
temperature, T 50 o, which is not equivalent to the gas tem- 
perature at r 5 oo for non-isothermal mass distributions [for 
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Figure 1. Projected mass- weighted entropy (top), surface gas mass density (middle), and projected mass- weighted temperature (bottom) 
maps of a randomly-selected galaxy group at z = in three different runs NoCool (left), ZCooLSFSN (middle), and ZCooLSF_SN_AGN 
(right). Each map is 10 h~ 1 Mpc on a side and centered on the group, which has a total mass of M500 ~ 3.3 X 10 13 Mq. The solid white 
circle delineates rsoo- Feedback by AGN has ejected high entropy gas from the high-z progenitors of the group. 



an NFW distribution, T(rsoo) ~ O.6T500, assuming hydro- 
static equilibrium and that gas traces dark matter]. There- 
fore, 5500 depends only on the halo mass, which is dominated 
by dark matter, and is insensitive to the thermodynamic 
state of the gas. 

For reference we show in Fig. 2 the 16th and 84th 
percentiles of the Chandra entropy profiles of Sun et al. 
(2009) [approximately 20 groups with a median mass of 
log 10 (M5oo/Mo) ss 13.5; note we have used here only the 



20 lowest mass groups of Sun et al. so as to achieve a similar 
median mass to that of the simulated groups]. Also shown 
(dotted line) is the power-law fit of VKB05 to the entropy 
profiles of a sample of groups and clusters simulated with 
non-radiative physics (see below). 

Since we are interested primarily in investigating what 
sets the global thermodynamic properties of the IGrM, we 
will focus our attention on quantities such as the median gas 
particle entropy of groups. Typically, the median entropy 
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ZCool_SF_SN_AGN 
VKB05 



0.1 1.0 

r / r 500 

Figure 2. Median 3D mass-weighted entropy profiles for a se- 
lection of OWLS runs for groups with masses between 13.25 ^ 
log lo (M 5 oo/M ) ^ 14.25. All other OWLS runs not shown (see 
Schaye et al. 2010) fall in between NoCool and ZCooLSF, in 
terms of their scaled entropy at intermediate radii. For reference, 
the hatched region represents the observational data of Sun et 
al. (2009) and the dotted line is the power-law fit of VKB05 to 
r > 0.3r5oo for a suite of non-radiative cluster simulations. Note 
that at intermediate radii the runs without SN feedback (i.e., 
ZCooLSF and PrimCooLSF) yield higher entropy than the runs 
that include SN feedback, while the run that includes AGN feed- 
back (ZCooLSF_SN_AGN) yields a similar entropy as the run 
with metal-line cooling and without feedback (ZCooLSF). 



corresponds to the entropy at intermediate radii of ~ 0.4 — 
0.6r5oo- From Fig. 2 we see that over this range in radii the 
simulations show a wide range of entropy values, spanning 
nearly an order of magnitude in S/ Ssoo ■ Of course, not all of 
the models we have included are physically reasonable, but 
several of them represent limiting cases which can be used to 
isolate the importance of certain subgrid physical processes, 
and the large differences between the models make it easier 
to do so. 

For example, NoCool, which is effectively a non- 
radiative run (note the agreement between the NoCool run 
and the power-law fit of VKB05, who fit approximately to 
the range 0.3 — 1.5rsoo), yields the lowest entropy of any of 
the runs. This is not unexpected, since the only process by 
which the entropy may be raised in this model is through 
gravitational shock heating. Since gravity has no preferred 
scale, the resulting shape and normalisation of the entropy 
profile does not depend on system mass, which is why it is 
sometimes referred to as the 'self-similar' profile. Note that 
within ~ 0.2r5oo the entropy profile deviates from a pure 
power-law behaviour, flattening into a core. This was also 
found by VKB05 (see also Frenk et al. 1999) for both adap- 
tive mesh refinement (AMR) and SPH non-radiative cosmo- 
logical simulations. The physical origin of the core appears 
to be related to the efficiency of energy exchange between 



gas and dark matter (McCarthy et al. 2007; see also Mitchell 
et al. 2009). 

ZCooLSF represents another interesting limiting case, 
as it includes cooling via both thermal bremsstrahlung and 
(metal) lines but no feedback. This model yields very high 
entropies at intermediate radii. That simulations without 
feedback can yield entropies in excess of the self-similar 
model has been noted previously (e.g., Muanwong et al. 
2002; Dave et al. 2002) and explained, at least qualitatively, 
by Knight & Ponman (1997) and by Bryan (2000), Voit & 
Bryan (2001), and Voit et al. (2002) (hereafter collectively 
referred to as VB). These authors argued that radiative cool- 
ing acts to selectively remove the lowest entropy gas from the 
hot phase, turning it into cold gas and stars. In the language 
of VB, radiative cooling 'truncates' the entropy distribution 
of the gas (Voit & Bryan 2001 and Voit et al. 2002 also noted 
the possibility that feedback could truncate the entropy dis- 
tribution, by removing the lowest entropy gas). As a result 
of this truncation, the mean entropy of the remaining hot 
gas increases with tim«0- This hot gas, which was initially 
at large radii, due to its buoyancy, would then flow inwards, 
approximately adiabatically, to fill the pressure decrement 
left by the gas that was converted into cold (pressure-less) 
form. As Fig. 2 emphasises, this results in a large amount 
of excess entropy in the remaining gas. We will show below 
that this is because of the unrealistically large fraction of 
baryons that are converted into stars in this run. 

The results of PrimCooLSF, which is identical to 
ZCooLSF except that the cooling rates are calculated as- 
suming primordial abundances, are also consistent with the 
VB truncation hypothesis. This run also yields entropy in 
excess of the NoCool (self-similar) result, but not to the same 
degree as ZCooLSF. This can be understood in terms of the 
reduced efficiency of cooling due to the neglect of metal-line 
emission in PrimCooLSF. The reduced cooling rate results 
in a smaller truncation of the entropy distribution of the hot 
gas which, in turn, leads to a smaller rise in the entropy. 

The relatively large difference in the resulting entropy 
distributions of ZCooLSF and PrimCooLSF highlights the 
importance of metal-line cooling in the simulations. Mod- 
els without metal-line cooling will seriously underestimate 
the fraction of baryons in stars relative to models that do 
include metal-line cooling. Comparisons of such models to 
observations can therefore quite easily lead to spurious con- 
clusions about the required efficiency of feedback processes. 
Note also that the difference in the resulting entropy distri- 
butions between runs with and without metal-line cooling 
is even larger when feedback is included (compare Prim- 
CooLSF _SN and ZCooLSFSN), since more particles are en- 
riched due to expulsion of metals from star-forming regions 
by SNe, particularly at high redshift. Unlike the OWLS runs, 
most simulations presented in the literature calculate cool- 
ing rates assuming primordial abundances only or, if they 
do include metals, they do so by assuming abundances that 
are a fixed fraction of solar. 

The other runs plotted in Fig. 2 (ZC00LSF.SN, Pnm- 



2 The entropy of an individual parcel of gas does not increase as 
a result of cooling (it decreases), but the mean/median entropy 
of the overall hot phase increases as a result of the removal of the 
lowest entropy gas from the hot phase over time due to cooling. 
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Figure 3. Gas mass (top), stellar mass (middle), and total baryon 
(bottom) mass fractions within r^oo as a function of group mass, 
M50O; for the OWLS runs investigated in this study. The mass 
fractions have been normalised to the universal baryon fraction, 
/i, = tlb/Qm. The dashed line in the bottom panel represents the 
median baryon fraction of the high mass systems in the SPH 
non-radiative Millennium Gas Simulation (Crain et al. 2007), 
while the dotted lines represent the 10th and 90th percentiles 
from that study. Variations in the included subgrid physics can 
have a large effect on the predicted stellar and gas mass frac- 
tions but the total baryonic fraction is always near the universal 
value. The one exception is the run that includes AGN feedback 
(ZCooLSF_SN_AGN), for which the total baryon fraction is a 
steep function of halo mass. Energy injection from BHs is able to 
eject gas from groups, particularly at high redshift. 



CooLSFSN, ZC00LSF.SN.AGN) also include feedback pro- 
cesses, so determining the relative importance of cooling 
and feedback in producing the excess entropy is not pos- 
sible based on Fig. 2 alone. Furthermore, we have not yet 
explicitly demonstrated that it is indeed the VB trunca- 
tion mechanism that is responsible for the excess entropy 
of the hot gas in the models without feedback (ZC00LSF 
and PrimCooLSF). In the remainder of this section we will 
exploit the Lagrangian nature of the OWLS simulations 
to deduce which processes are most important for raising 
the entropy of the IGrM in the simulations. Of particu- 
lar interest will be isolating the effects the AGN feedback, 
since only ZCooLSF_SN_A GN reproduces the observed low 
baryon fractions of galaxy groups (see Paper I). 



3.2 Excess entropy in runs without SN and AGN 
feedback 

Before discussing the Lagrangian entropy distributions and 
histories, we present Fig. 3, which shows the gas, stellar, 
and total baryon mass fractions within r^oo as a function 
of mass, M500, of the simulated groups in all of the runs 
we consider in the present study. This figure is useful for 
the interpretation of the results presented below and will be 
referred to frequently. 



3.2.1 Present-day entropy distributions 

In Fig. 4 we plot the z — sorted gas entropy as a function of 
the included gas mass fraction for the NoCool, PrimCooLSF 
and ZC00LSF runs. The sorted entropy is calculated for each 
group by ordering the hot (T > 10 5 K) gas particles within 
fsoo by their entropies and calculating the total mass of gas, 
M gas (< S), with entropy lower than S. The motivation for 
examining S(M gas ), rather than S(r), is that cooling modi- 
fies S(M glls ) in a very simple way: it selectively removes the 
gas with entropy lower than the cooling threshold (see VB; 
see also below). The resulting radial entropy distribution, 
by contrast, depends not only on the resulting entropy, but 
also on the depth of the gravitational potential well, since 
the gas will tend to re-establish hydrostatic and convective 
equilibrium. For this reason, it is conceptually much easier 
to deduce the effects of cooling on the IGrM using the sorted 
gas entropy as a function of the included gas mass fraction. 
The solid curves represent the median relations for the three 
runs for groups with 13.25 ^ log lo (Msoo/Af0) ^ 14.25. 

Note that for NoCool the gas mass fraction within rsoo 
is approximately 0.9 times the universal baryon fraction, fb 
(the blue curve ends at M gas (< S) « 0.9ft, Msoo). This is in 
accordance with the results of previous non-radiative SPH 
simulations which calculated the baryon fraction within a 
similar fixed aperture (e.g., Eke et al. 1998; Frenk et al. 
1999; Kravtsov et al. 2005; Crain et al. 2007). The gas mass 
fractions for ZC00LSF and PrimCooLSF, however, are much 
lower - approximately 0.3 and 0.4 times the universal baryon 
fraction, respectively (see top panel of Fig. 3). In these runs 
most of the baryons have been able to cool significantly 
and form stars (middle panel of Fig. 3). However, the to- 
tal baryon fraction (stars + gas) is very similar to that of 
NoCool, as can be seen from the bottom panel of Fig. 3. 
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Figure 4. The sorted gas entropy as a function of included gas 
mass fraction for the NoCool, ZCooLSF and PrimCooLSF runs 
(solid blue, green, and red curves, respectively). This is calculated 
for each group by ordering the gas particles within rsoo by their 
entropies and calculating the total mass of gas with entropy lower 
than S, A-f gaa (< S). The solid curves represent the median rela- 
tions for the 3 runs for groups with 13.25 ^ log 1() (M5oo/A^©) ^ 
14.25. The dashed green and red curves represent the sorted en- 
tropy as a function of enclosed gas mass fraction for the No- 
Cool run when all gas with entropy lower than that which en- 
closes / gas = /*,zCool_SF and /.,PrimCoolJ3F> respectively, is re- 
moved ("truncated"). This operation yields the expected entropy 
distribution after cooling has selectively removed the lowest en- 
tropy/shortest cooling time gas. Overall, the truncated entropy 
distributions are similar to those obtained directly from ZCooLSF 
and PrimCooLSF, indicating that the excess entropy in these runs 
is indeed largely the result of selective removal of low entropy gas, 
as proposed by VB. See text for further discussion. 



Kravtsov et al. (2005) find very similar trends for their sim- 
ulations that include radiative cooling and star formation 
(and inefficient feedback). 

A comparison of the solid curves reveals excess entropy 
at fixed gas mass fraction in the runs with cooling rela- 
tive to the baseline NoCool run. This plot is similar to the 
trend seen in Fig. 2 because the buoyancy of the gas re- 
sults in a monotonic mapping between r and M gas (< S). 
According to the model proposed by VB, cooling selectively 
removes the lowest entropy gas from the hot phase, which 
results in a shift in the S(M gaa ) curves to the left in Fig. 
4. To see whether this idea works in detail, we truncate the 
NoCool S(M g!lB ) distribution (i.e., shift to the left) using 
the calculated stellar fractions from the simulated groups in 
ZCooLSF and PrimCooLSF (dashed and dot-dashed curves, 
respectively). This is done by removing all gas with entropies 
lower than the entropy that encloses / gas = /» . The remain- 
ing gas mass fraction is then just the original gas mass frac- 
tion minus the stellar mass fraction. This procedure yields 



the expected S(M gas ) distribution after cooling has removed 
the lowest entropy gas from the hot phasfl 

Overall, the entropy distributions derived by truncat- 
ing the NoCool distribution using the stellar fractions from 
ZCooLSF and PrimCooLSF are similar to the actual entropy 
distributions of the hot gas obtained from these runs. This 
indicates that truncation due to cooling is the likely explana- 
tion for the bulk of the excess entropy in these runs. But note 
that a shortfall of a factor of approximately 1.5 at interme- 
diate gas masses is present, implying that the VB truncation 
mechanism is not a perfect description of what happens in 
the simulations. A plausible explanation for this shortfall is 
that the remaining hot gas in ZCooLSF and PrimCooLSF is 
(gravitationally) shocked to a higher degree than the same 
gas in the NoCool run because it is able to fall further into 
the gravitational potential of the group, and therefore shock 
at a higher Mach number, due to the reduced baryon frac- 
tion (and thermal pressure) of the group (e.g., McCarthy 
et al. 2007). Indeed, as we discuss below, the gas that con- 
stitutes the IGrM in ZCooLSF and PrimCooLSF is largely 
located at r > rsoo in NoCool, implying thermal pressure 
has held the gas up in the latter. 



3.2.2 Lagrangian thermal histories 

To test these ideas, we calculate for each group the median 
entropy, density, and temperature as a function of redshift 
for all the gas that ends up hot (T > 10 K) and within rsoo 
at the present day (i.e, that which constitutes the IGrM). 
In Fig. 5 we show the median histories for galaxy groups 
in the mass range 13.25 ^ log 10 (M 50 o/M©) ^ 14.25. The 
solid blue and red curves represent NoCool and ZCooLSF, 
respectively. 

Focusing first on the entropy histories (top panel), we 
see that prior to z ~ 2 the two runs are quite similar; they 
show a gradual rise in the gas entropy due to both gravi- 
tational shock heating (i.e., particles falling into potential 
wells, which are becoming deeper with time) and UV pho- 
toheating. The sharp rise in entropy at z — 9 corresponds 
to reionisation, which is assumed to occur everywhere in 
the simulation volume simultaneously. At late times, the en- 
tropy histories diverge, such that by the present day (z = 0) 
there is approximately a factor of 4 difference in the medi- 
ans, consistent with the z — entropy profiles plotted in 
Fig. 2. However, it is important to note that by selecting 
the hot gas within rsoo at z = for the two simulations, 
we are not selecting identical sets of particles. The gas that 
ends up constituting the IGrM in one simulation may differ 
from that which ends up constituting the IGrM in another 
simulation. This implies that we are comparing the ther- 
mal histories of different gas 'parcels'. However, by using 
the unique IDs of the gas particles, we can select identi- 
cal sets of particles in the simulations and compare their 



3 More precisely, this represents an upper limit to the entropy 
distribution after cooling has selectively removed the lowest en- 
tropy gas, as it has been implicitly assumed that any gas that 
has not cooled out of the hot phase flows adiabatically inward to 
replace the gas that has cooled out. In reality, however, some of 
the remaining hot gas would have its entropy reduced by cooling 
as its density increases. 
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Lagrangian thermal histories (note by using massless tracer 
particles, similar analyses can also be applied to mesh-based 
simulations, e.g., Vazza 2010). 

The dashed blue curve in the top panel of Fig. 5 
shows the median entropy history of gas in the NoCool 
(non-radiative) run that ends up constituting the IGrM in 
ZCooLSF (i.e., is hot and within rsoo at z = in ZCooLSF). 
(But note that there is no requirement that this gas ends up 
as part of the IGrM in the NoCool run itself.) By comparing 
this curve to the solid red curve, we are making a comparison 
of identical sets of particles in two runs. As can clearly be 
seen, the histories are remarkably similar, with the median 
entropies differing only by a factor of 1.5 or so at the present 
day. This agrees with the difference between the solid and 
dashed green curves seen in Fig. 4. Since we selected and fol- 
lowed back in time the same set of particles, which are those 
particles that end up hot and in groups at z — in ZCooLSF, 
the similarity of the histories firmly establishes that the VB 
truncation mechanism is responsible for the bulk of the ex- 
cess entropy in ZCooLSF. We can make this statement as we 
have selected only those particles that survived truncation. 
The similarity of the histories implies that the reason for the 
large difference in the z = (non-truncated) entropy distri- 
butions is a result of cooling selectively-removing the lowest 
entropy gas. As also predicted by VB, the remaining hot 
gas flows inward to replace the truncated low entropy gas. 
The median radius of the hot gas within rsoo of ZCooLSF is 
approximately 0.7r5oo, whereas the same set of gas particles 
have a median radius of approximately 1.3r5oo in NoCool. 
Hot gas in NoCool is unable to cool and its pressure prevents 
gas outside of 7-500 from falling inward. 

We can also perform the inverse test, which is to select 
from the NoCool run those particles that end up hot and 
within rsoo of galaxy groups and show their history in the 
ZCooLSF run. This is represented by the dashed red curve 
in the top panel of Fig. 5. Of course, a large fraction of the 
gas that ends up hot and within rsoo in the NoCool run will 
be able to cool and form stars in the ZCooLSF run. The 
dashed red curve therefore represents the median entropy of 
the gas which does not get converted into stars (therefore 
the number of gas particles we use to calculate the median 
entropy changes with time for this curve). The dot-dashed 
green curve represents the fraction of gas particles (that end 
up hot and within rsoo in the No Cool run) that are converted 
into star particles in the ZCooLSF run. 

At early times, z > 5, cooling is unimportant for the gas 
in the ZCooLSF run that ends up constituting the IGrM in 
the NoCool run and, as a result, the entropy history of this 
population of particles in the ZCooLSF run is similar to that 
of the same gas in the NoCool run. At 2 < 5, cooling becomes 
rapidly more important and we see a departure between the 
populations in the two runs, with the gas in the ZCooLSF 
run having stalled its entropy production. Not coinciden- 
tally, this departure corresponds closely to the time when 
the gas begins to be rapidly converted into stars. Eventu- 
ally, by 2 ~ 1 all of the low entropy gas is converted into 
stars (see dot-dashed green curve). The entropy history of 
the remaining hot gas then becomes very similar to that rep- 
resented by the solid red curve, since these are now almost 
identical sets of particles. 

All of the above information was gleaned from an analy- 
sis of the entropy history alone, but it is worthwhile to briefly 
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Figure 5. Top: The solid blue and red curves represent for No- 
Cool and ZCooLSF, respectively, the median entropy history of 
gas that ends up hot and within rsoo at the present day for 
groups in the mass range 13.25 ^ log^A/soo/A^©) ^ 14.25. The 
dashed blue curve represents the median entropy history of gas 
in NoCool that ends up hot and within rsoo in ZCooLSF. The 
dashed red curve represents the median entropy history of gas in 
the ZCooLSF run that ends up hot and within rsoo in the No- 
Cool run. The dot dashed green curve represents the fraction of 
group gas particles within rsoo in NoCool that have been con- 
verted into star particles in ZCooLSF. Middle: The analogous 
plot for the gas density. The gas density has been normalised by 
the 'virial gas density' psoo( z = 0) = 500/j,p C rit(2 = 0). Bottom: 
The analogous plot for the gas temperature. The temperature 
has been normalised by the 'virial temperature' Tsoo(z = 0) = 
GM5ooM m ff /(2r5oofcs). When identical sets of particles are se- 
lected in the two runs (compare the solid red and dashed blue 
curves) the histories are similar. This implies that the excess en- 
tropy in the ZCooLSF run is mostly due to the selective removal 
of the lowest entropy gas via cooling and star formation. 
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consider the density and temperature histories as well. In the 
middle and bottom panels of Fig. 5 we see that gas that is 
destined to constitute the IGrM in the ZCooLSF run (solid 
red) is typically cooler (at least until recently) and less dense 
than that which will constitute the IGrM in the NoCool run 
(solid blue). Given the similarity of the entropy histories of 
the two runs until z ~ 2, this suggests the IGrM in the No- 
Cool run was, on average, compressed (adiabatically) more 
strongly at early times, presumably because it was typically 
located deeper within the potential wells of small haloes. Gas 
that is destined to constitute the IGrM in the ZCooLSF run, 
by contrast, is that which preferentially avoided significant 
cooling due to its lower density, presumably because it was 
located more on the periphery of overdensities at early times. 
This is confirmed by examining the dashed blue curve. 

It is also worth pointing out that, had we focused on 
the density and/or temperature alone without consideration 
for the entropy, we might have reasonably (but wrongly) 
deduced that cooling was directly important at very early 
times (z >5), since the densities and temperatures of the gas 
in the two runs differ at that time. However, the similarity of 
the entropies of the gas demonstrates that this is due purely 
to differences in the (adiabatic) compression of the gas at 
early times. 

We have also examined the thermal history of the other 
run without feedback (PrimCooLSF) and find that it is very 
similar to that of ZCooLSF and confirms the above conclu- 



3.3 Excess entropy in runs with SN feedback 

In Section 3.2 we demonstrated that in runs without feed- 
back the excess entropy at intermediate radii is largely a 
consequence of the selective removal of low entropy gas due 
to cooling and star formation. Runs that include SN feed- 
back also show evidence for excess entropy relative to the 
self-similar model (see Fig. 2), but it is not immediately ob- 
vious in this case whether it is cooling/star formation or SN 
heating (or both) that is responsible for the extra entropy. 
A comparison of the thermal history of the groups in runs 
with SN feedback with that of NoCool can resolve this issue. 

In Fig. 6 we show the median entropy, density, and 
temperature histories for galaxy groups in the mass range 
13.25 s; log 10 (M 50 o/M Q ) 14.25. The solid blue and red 
curves represent NoCool and ZCool.SF.SN, respectively. Fo- 
cusing first on the entropy (top panel), at the present day 
(z = 0) there is approximately a factor of 2 difference be- 
tween the medians in the two runs. This is smaller than the 
difference we found between NoCool and ZCooLSF (Fig. 5). 
This is already a clue that the addition of SN feedback does 
not significantly raise the entropy of the IGrM. 

Selecting from NoCool only those particles which end 
up being hot and within r^oo in ZCool.SF.SN (dashed blue 
curve), we see a remarkably similar thermal history to this 
same particle population in ZCool.SF.SN (solid red curve). 
This immediately tells us that SN feedback has not sig- 
nificantly raised the entropy of gas that ends up hot and 
within groups in ZCool.SF.SN. It is conceivable, though, 
that some of the gas could have been heated significantly 
by SN feedback and, as a result, did not end up within rsoo 
in ZCool.SF.SN. However, the similarity of the dashed red 
curve (i.e., those gas particles in ZCool.SF.SN that end up 
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Figure 6. Same as Fig. 5 but for ZCool.SF.SN. The similarity of 
the solid red and dashed blue curves indicates that SN feedback 
has not significantly raised the entropy of the gas that consti- 
tutes the present-day IGrM in ZCool.SF.SN. The similarity of 
the dashed red and solid red curves indicates that SN feedback 
has not strongly heated any gas that may have been ejected ei- 
ther. 



hot and within groups in NoCool) and the solid red curve 
(gas that ends up hot and in groups in the ZCool.SF.SN 
run) demonstrates that SN feedback did not significantly 
raise the entropy of any gas ever associated with the groups 
in ZCool.SF.SN. 

As in the case of the runs without SN feedback, we 
thus find that the primary cause of excess entropy in 
ZCool.SF.SN is the selective removal of low entropy gas. As 
already noted, the amount of excess entropy in this run is 
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smaller than for runs without feedback. In addition, the stel- 
lar fraction in this run, /*(rsoo) ~ 0.4/b, is smaller than that 
for the runs without SN feedback (ss 0.6/f, for ZCooLSF; see 
Fig. 3). The role that SN feedback plays, therefore, is to 
prevent some of the low entropy gas from turning into cold 
gas/stars which implies that the entropy distribution is not 
truncated to the same degree as in the runs with no effec- 
tive feedback. Consequently, the entropy of the IGrM is not 
raised by the same amount. And feedback from SN-driven 
winds does not significantly heat the IGrM. It is worth not- 
ing here that since the SN feedback is implemented in kinetic 
form in the OWLS simulations (see Dalla Vecchia & Schaye 
2008), it does not suffer from the standard problem asso- 
ciated with most thermal implementations of SN feedback; 
i.e., that the gas radiates away its gained energy almost im- 
mediately. 

The density and temperature histories yield a similar 
story to what we found by analysing the runs with no feed- 
back; i.e., the gas that constitutes the IGrM at the present- 
day in the ZCool.SF.SN run is that which escaped significant 
cooling because it was of lower density and temperature at 
early times. Again, this is presumably because it was located 
in the outskirts of haloes at early times. 

The thermal history of the other run with SN feedback 
we consider (i.e., PrimCool.SF.SN) is very similar to that 
of ZCool.SF.SN, in the sense that SN feedback does not 
noticeably raise the entropy of the IGrM. Instead, it acts 
primarily to prevent some of the lowest entropy gas from 
turning into stars/cold gas. However, the lack of metal-line 
cooling in this run greatly reduces the fraction of baryons 
converted into stars, resulting in a smaller truncation of the 
entropy distribution and entropy profiles that clearly are at 
odds with observations (see Fig. 2). 

3.4 Excess entropy in a run with AGN feedback 

The OWLS runs we have investigated so far, particularly 
those which include the important process of metal-line 
cooling, produce galaxy group stellar mass fractions that 
are much larger than observed. In particular, current esti- 
mates suggest that /*(rsoo) is between 0.1-0.2 for nearby 
groups (Lin & Mohr 2004; Balogh et al. 2008; McGee & 
Balogh 2010), depending on what fraction of the total light 
is assumed to be in a diffuse component (see, e.g., Gonza- 
lez et al. 2007; McGee & Balogh 2010). This implies that 
ZCool.SF.SN and ZCooLSF have too high stellar mass frac- 
tions by factors of « 3 and 4, respectively (see Fig. 3). This 
is the so-called 'cooling crisis' of cosmological hydrodynamic 
simulations that invoke inefficient feedback (e.g., Balogh et 
al. 2001). 

Based on the above results, any simulation that has too 
high a stellar mass fraction will overestimate the degree to 
which the entropy of the IGrM will be raised by cooling/star 
formation, because it overestimates the degree to which the 
self-similar entropy distribution has been truncated. It is 
interesting to note that of the runs we have examined so 
far, the only ones with excess entropy comparable to what 
is observed (ZCooLSF, PnmCooLSF, ZCool.SF.SN) signif- 
icantly violate the observed stellar mass fraction of groups. 
It is also interesting that while the model that includes SN 
feedback but calculates cooling rates assuming primordial 
abundances (PrimCool.SF.SN), predicts a stellar mass frac- 
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Figure 7. Same as Fig. 5 but for ZCooLSF.SN.AGN. The simi- 
larity of the solid red and dashed blue curves indicates that feed- 
back from BHs has not strongly heated the gas that ends up con- 
stituting the present-day IGrM. However, the deviation of the 
dashed red and solid red curves at z ~ 2 implies that AGN feed- 
back has significantly heated gas which did not end up in groups 
in ZCool.SF.SN.AGN but did in NoCool (i.e., ejected gas). 



tion only slightly larger than observed, the predicted entropy 
is much lower than observed^]. It therefore does not appear 
possible to simultaneously match the observed thermal prop- 



4 Of course PrimCool.SF.SN ignores the important effects of 
metal-line cooling, so it does not represent a physically reason- 
able model. We note, however, than many cosmological simula- 
tions found in the literature adopt primordial abundances when 
calculating cooling rates. 
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Figure 8. The fraction of gas removed (truncated) from the hot 
IGrM due to star formation in ZCooLSF or star formation plus 
ejection in ZCool.SF.SN.AGN as a function of group mass. Also 
shown is the fraction removed from the hot IGrM due to star 
formation alone in ZCool.SF.SN.AGN. The reason the two runs 
have similar present-day entropy at intermediate radii is that both 
target the same low-entropy gas for truncation. 



erties of the IGrM and the observed stellar mass fractions 
of groups in the context of a model where the entropy is 
"raised" (by truncation of low entropy gas) via cooling/star 
formation alone. This likely indicates that some form of 
highly efficient feedback (in terms of raising the entropy) 
is required in order to explain the excess entropy and stellar 
mass fractions of galaxy groups simultaneously. 

We demonstrated in Paper I that model 
ZCooLSF.SN.AGN, which includes AGN feedback, yields 
the correct stellar mass fraction and entropy distribution 
for the IGrM as a function of halo mass, in addition to 
matching a wide range of additional properties. What is the 
physical process that sets the entropy distribution of the 
IGrM in this model? Four possible scenarios are as follows: 

• BH feedback raises the entropy of the gas which, in 
turn, increases its cooling time and reduces its star forma- 
tion efficiency. Since the peak of the star formation (and 
presumably BH formation) history of the universe occurs 
at z m 2 — 3, and present-day groups contain most of the 
galaxies, this scenario would imply that the excess entropy 
was generated at high redshift, well before the formation 
of the IGrM. This is another way of saying the IGrM was 
'preheated' (e.g., Kaiser 1991; Evrard & Henry 1991). 

• A second possibility is that the entropy was raised much 
more recently through direct heating by BH^. 

• A third possibility is that BH feedback targets only 



5 Although McCarthy et al. (2008) argued that such internal 
heating is energetically prohibitive due to the high density of the 
gas after the system has collapsed. 



the lowest entropy /shortest cooling time gas for feedback, 
ejecting it from the low-mass progenitors of present-day 
groups, while leaving the long cooling time gas unaltered 
(i.e., the long cooling time gas ends up forming the present- 
day IGrM). Here we do not distinguish whether the low en- 
tropy gas was first heated and then buoyantly floated out 
of the system or if it was launched balistically (i.e., given a 
large velocity kick) from the system. But note the entropy 
gain/velocity kick must be sufficiently large so that the gas is 
not re-accreted as the group halo grows. Some tuning of this 
mechanism is needed, as the gas is recaptured by massive 
clusters (i.e., their baryon fractions are near the universal 
value) . 

• Low entropy gas is ejected from present-day groups, as 
proposed recently by B08. To be more accurate, B08 actually 
allowed for ejection from the low mass-progenitors of groups 
as well. However, this gas is re-accreted by the group at 
late times if the energy imparted to the gas is less than 
the binding energy of the present-day group. Therefore, the 
energetics required to eject the gas are set by the depth of 
the potential well of the group at z — 0, which is why we 
refer to B08 as a 'present-day ejection' model. 

The last two scenarios differ from the first two in that 
the excess entropy in groups is not achieved by heating of the 
IGrM but instead by the removal of the lowest entropy gas, 
as proposed by VB, modulo that the truncation results from 
the combination of cooling/star formation and ejection. 

An examination of the thermal history of 
ZCool.SF.SN. AGN can distinguish between these four 
proposed scenarios (i.e., early heating, late heating, 
early ejection, and late ejection, respectively). In Fig. 
7 we show the median entropy, density, and temper- 
ature histories for galaxy groups in the mass range 
13.25 < log 10 (M 50 o/M Q ) < 14.25. The solid blue and 
red curves represent NoCool and ZCool.SF.SN.AGN, 
respectively. Focusing first on the entropy (top panel), at 
the present day (z — 0) there is approximately a factor 
of 4 difference between the medians in the two runs, 
which is similar to the difference between ZCooLSF and 
NoCool (Fig. 5). We discuss below whether or not this is a 
coincidence. 

Selecting in the NoCool run only those particles that 
end up hot and within rsoo in ZCool.SF.SN.AGN (dashed 
blue curve), we see a similar thermal history as this same 
particle population has in ZCool.SF.SN.AGN (solid red 
curve). This is somewhat surprising as it indicates that feed- 
back from supermassive BHs has not substantially raised the 
entropy of gas that ends up constituting the IGrM in the run 
with AGN feedback This immediately rules out the first 
two scenarios we outlined above, which are based on the idea 
that the excess entropy in groups is a result of heating of 
the gas that ends up in the IGrM. 

Of course, some amount of gas could have been heated 
significantly by AGN feedback and, as a result, ended up 
outside of rsoo in ZCool.SF.SN.AGN. Indeed, the reduced 



6 Note a slight offset is present between the solid red and dashed 
blue curves. However, as we discuss in Section 4.2, this is most 
likely not a consequence of AGN feedback but, rather, amplified 
gravitational shock heating due to the reduced thermal pressure 
of the IGrM in the groups in the run with AGN. 



How quasar winds solve over-cooling in groups 13 



baryon fractions of groups in the ZCooLSF_SN_AGN run 
relative to the universal value (see Fig. 3) already indicates 
that this is the case. This is reflected in the entropy history 
as a departure of the dashed red curve (i.e., those particles 
in ZCooLSF_SN_AGN that end up hot and within groups in 
NoCool) from the solid red curve (gas that ends up hot and 
in groups in ZCooLSF.SN.AGN) at z w 2. BH feedback has 
therefore raised the entropy of some of the gas that would 
have ended up in groups if not for this boost. This is contrary 
to what we found for the runs with only SN feedback (see 
Fig. 6). 

Thus, the way in which the entropy of the IGrM is raised 
by feedback from BHs is not by heating of gas destined to 
end up in groups, but instead by ejecting low entropy gas at 
early times. We noted above that the degree to which the 
entropy is raised in ZCooLSF_SN_AGN is similar to that in 
the run that includes metal-line cooling and star formation 
but no significant sources of feedback (ZCooLSF). At first 
sight this would appear to be a coincidence. However, this 
is exactly what one expects if the excess entropy is due pri- 
marily to the removal of just the lowest entropy gas (which 
happens via star formation in ZCooLSF and star formation 
plus ejection in ZCooLSF.SN.AGN). We therefore expect 
that the sum of the stellar mass and ejected gas mass frac- 
tions for ZCooLSF_SN_A GN to be similar to the stellar mass 
fraction in ZCooLSF, which is something that can be explic- 
itly checked in the simulations. 

In Fig. 8 we show the fraction of gas that has been 
truncated as a function of M500 for ZCooLSF_SN_AGN and 
ZCooLSF. The truncation fraction in ZCooLSF is simply 
the stellar mass fraction, whereas for ZCooLSF_SN_AGN it 
is the summation of the stellar mass fraction and the ejected 
gas mass fraction. We compute the ejected gas mass fraction 
for each group in the ZCooLSF_SN_AGN as the difference 
in the baryon fraction of the same group between ZCooLSF 
and ZC00LSF.SN.AGN. 

We see that over a large range in halo mass (M500 < 
5 x 1O 13 M0) the degree of truncation in the two simulations 
is indeed quite similar, providing strong evidence that it is 
indeed the lowest entropy gas that has been targeted by BH 
feedback. Any low entropy gas that escapes BH feedback 
ends up forming stars, but clearly it is ejection, rather than 
star formation, that dominates the truncation fraction, as 
can be quickly deduced by comparison to the stellar mass 
fractions of the ZC00LSF.SN.AGN run. 

Interestingly, at high halo masses the truncation frac- 
tions of the two runs begin to depart, with / tranc approach- 
ing /* for ZC00LSFJSN-AGN. This is not unexpected, as it 
becomes increasingly more difficult for the BH feedback to 
eject gas beyond rsoo, due to the increased binding energies 
of more massive systems (see also Fig. 3, which shows the 
total baryon fraction is approaching the universal value for 
these very massive systems). 



4 THE PHYSICAL NATURE OF AGN 
FEEDBACK 

We demonstrated in Section 3.4 that the way in which feed- 
back from AGN raises the entropy of the IGrM is through 
the ejection of low entropy gas, rather than heating of the 
gas that constitutes the present-day IGrM. This interesting 
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Figure 9. The timing of entropy injection due to AGN feedback 
in the ZCooLSF JSN^AGN run. Shown are the median 3D mass- 
weighted entropy profiles at z = for the AGN feedback model 
when the AGN have been switched off in the simulation at red- 
shift z = z g (see text). The median is calculated for groups in 
the mass range 13.25 ^ log 10 (JWsoo / ^q) C 14.25. For reference, 
the hatched region represents the observational data of Sun et 
al. (2009). Feedback from AGN-driven outflows has significantly 
heated much of the gas that had the lowest entropy at high red- 
shift, ejecting it from groups primarily at z > 1.5. 

result warrants further investigation. In particular, we would 
like to know how, when, and to what degree the ejected gas 
is heated by feedback from BHs. We address these questions 
below. We also briefly discuss the potential sensitivity of our 
conclusions to the adopted parameters of the AGN feedback 
model. 



4.1 When was the gas ejected/heated? 

In Fig. 9 we show the median radial entropy distribution 
of the gas for runs using the AGN model where AGN feed- 
back has been explicitly shut off in the simulation at redshift 
z — ZoB (i.e., we have re-run the AGN model a number of 
times, each time explicitly switching off AGN feedback at 
a different redshift). Switching off the AGN in the simu- 
lation provides a simple tool for deducing the importance 
of AGN feedback subsequent to z ff. The simulation with 
metal-dependent cooling, star formation, and feedback from 
SNe but not from AGN (ZC00LSF.SN) would correspond to 
the case of z g — 00 in this terminology. 

Comparing the dashed black curve (z ff = 00) and solid 
blue (z ff = 1.5) and red (AGN never shut off) curves at 
intermediate radii, we see that the bulk of the excess entropy 
relative to the self-similar profile is produced at z > 1.5, 
since the solid blue curve is about half way between the 
dashed black and solid red curves. Note that this is likely a 
lower limit for the redshift where AGN feedback raises the 
entropy of the IGrM, since excess entropy can be produced 
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Figure 10. Gas mass (top), stellar mass (middle), and total 
baryon (bottom) mass fractions within rsoo as a function of group 
mass, M500, for the AGN feedback model when the AGN have 
been switched off in the simulation at redshift z = z g . The mass 
fractions have been normalised to the universal baryon fraction, 
ft = Q(,/n m . Feedback from AGN-driven outflows has signifi- 
cantly heated much of the lowest entropy gas, ejecting it from 
groups primarily at z > 1.5. 



even after the AGN is switched off via truncation due to 
radiative cooling (as in Section 3.3). 

We show the effect of switching off the AGN on the 
gas, stellar, and baryon mass fractions within r 5 oo at z = 
in Fig. fO. Consistent with the entropy profiles shown in 
Fig. 9, we see here that the bulk of the ejection of gas (and 
the prevention of star formation) through AGN feedback 
occurred at z > f.5, as the blue filled circles (represent- 
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Figure 11. The median entropy history of ejected gas in 
ZCooLSF_SN_AGN. (An ejected gas particle is one that is lo- 
cated within rsoo in NoCool but is outside this radius in 
ZC00LSFSN-AGN.) Feedback from AGN-driven outflows has 
significantly heated much of the gas that had the lowest entropy 
at high redshift, ejecting it from groups primarily at z > 1.5. 



ing the case z a a = 1.5) are closer to the solid filled circles 
(ZC00LSFSN-AGN) than to the black crosses (z off = 00) 
in all three panels shown in Fig. 10. The AGN model solves 
the overcooling problem common to many previous cosmo- 
logical simulations because the feedback is efficient at high 
redshift, which is when the bulk of the star formation occurs 
(e.g., Nagai & Kravtsov 2004). 

A more rigorous, though perhaps less intuitive, determi- 
nation of the epoch during which the gas was heated/ejected 
can be deduced from the median entropy history of ejected 
gas in ZCoolJSFJSN-AGN, which is shown in Fig. 11. Also 
shown (solid black) is the median entropy history of the gas 
that ends up constituting the IGrM in this run. 

At early times (z > 4) the entropy history of ejected gas 
is similar to that of gas which is not ejected. This indicates 
that feedback from BHs is unimportant prior to z ~ 4. How- 
ever, over a relatively small range of redshift (2 < z < 4), 
we see a sharp departure of the ejected gas from the gas des- 
tined to constitute the IGrM, indicating that AGN feedback 
is very important during this time. This behaviour is what 
we naively expect given the mass accretion rate history of 
the model (see Fig. 7 of Booth & Schaye 2009, which shows 
that the BHs gain most of their mass between z — 4 and 
z — 2). Observationally, this is the epoch when the quasar 
luminosity function peaks, and the energy injection rate is 
proportional to the luminosity. It is important to note that 
most of the BH growth occurs when the accretion rate is 
close to Eddington-limited, which is when the BH would 
be considered to be in 'quasar mode', rather than in 'radio 
mode'. That our simulations indicate that the entropy of the 
present-day IGrM is largely shaped by AGN 'quasar mode' 
feedback differs from previous studies (e.g., B08; Giodini et 
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al. 2010) and is one of our major results. Note that the time 
period over which the heating/ejection occurs also corre- 
sponds to the peak of the star formation rate density of the 
universe for this run (see Booth & Schaye 2009), indicating 
that there is a strong correlation between BH growth and 
star formation in the simulation, as is also observed (e.g., 
Boyle & Terlevich 1998; Kauffmann et al. 2003; Hfiring & 
Rix 2004). Indeed, this is presumably the reason why the 
simulation reproduces the z = black hole mass - stellar 
mass relation. 

The ordering of the curves in the right panel of Fig. 1 1 is 
also very interesting. It indicates that the lowest entropy gas 
is the first to be ejected and ultimately ends up gaining the 
most entropy. This is consistent with the expectation that 
the gas that is ejected early on is able to travel further away 
from the group by z = and will therefore have acquired 
more entropy through weak shocking as it interacts with the 
intergalactic medium. 

4.2 To what degree is the ejected gas heated? 

To more clearly demonstrate the effect of BH feedback on 
gas that did not end up in groups in ZCool.SF.SN.AGN h\it 
did end up in groups in NoCool (i.e, ejected gas), we present 
Fig. 12. For each group at z — we compute the distri- 
bution of the logarithm of the ratio of the particle entropy 
in ZCooLSF.SN.AGN to the particle entropy in NoCool for 
particles that end up either inside or outside of groups in 
ZCooLSF.SN.AGN. This is done on a particle- by-particle 
basis by exploiting the identical initial conditions and unique 
particle IDs of gas in the simulations. For comparison, wc 
also compute the analogous distributions for ZCool.SF. The 
curves in Fig. 12 represent the median histograms for groups 
in the mass range 13.25 < log lo (Msoo/M0) < 14.25. 

Fig. 12 shows that the z — distribution of en- 
tropy magnifications relative to NoCool for particles that 
end up inside groups in ZCool.SF.SN.AGN is similar to 
that of particles that end up either inside or outside of 
groups in ZCool.SF, which has no feedbaclfl This can only 
mean that the gas that ends up constituting the IGrM in 
ZCool.SF.SN.AGN has not been significantly heated by BH 
feedback, as discussed in Section 3.4. 

By contrast, the distribution of entropy magnifi- 
cations for particles that end up outside groups in 
ZCool.SF.SN.AGN (i.e., ejected gas) is shifted away from 
the other curves, implying that the gas has been strongly 
heated on average, and the distribution is skewed to higher 
ratios (i.e., a larger fraction of the particles have very high 
entropy ratios) . Typically, the entropy of an ejected gas par- 
ticle in the run with AGN feedback has had its entropy 
raised by a factor of ~ 10 over that in NoCool, although 
many ejected particles receive boosts of up to ~ 100 times 
(and higher). 



7 Note that the distributions for ZCool.SF are not centered on 
(in log space), as one might naively expect. As discussed in 
Section 3.2.1, this is most likely because the gravitational heating 
experienced by the IGrM in this run exceeds that of the same gas 
in NoCool due to the reduced baryon fraction of groups in the 
former, which allows gas to fall further into the potential wells of 
groups and therefore to shock to a higher entropy. 
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Figure 12. The distribution of the ratio of the particle entropy 
in ZCool.SF.SN.AGN (red curves) or ZCooLSF (blue curves) 
to the particle entropy in NoCool, both for particles that end 
up inside (dashed curves) and outside (solid curves) groups in 
ZCool.SF.SN.AGN or ZCooLSF. The similarity of the dashed 
red and dashed blue curves indicates that AGN feedback has not 
strongly heated the gas that ends up constituting the IGrM. How- 
ever, the solid red curve, which is shifted to higher entropy mag- 
nifications, indicates that the gas outside groups in the run with 
AGN feedback has been strongly heated and was ejected from the 
system as a result. 



4.3 How was the gas heated? 

Was the gas heated abruptly or gradually over time? Until 
now we have focused only on median entropy histories, which 
is the typical gas particle entropy as a function of time. The 
advantage of the median is that it gives an indication of what 
the population of particles as a whole is doing. The disad- 
vantage is that it effectively 'smooths' out the histories of 
individual particles. Thus, the median history is not partic- 
ularly useful for assessing if the heating process occurs grad- 
ually or through one abrupt episode. Gradual heating could 
happen, for example, if outflowing gas experiences numer- 
ous weak shocks from collisions with hydrostatic/inflowing 
gas in galaxies or with other outflowing gas (e.g., newly out- 
flowing gas could catch up with gas heated earlier), or if 
the gas re-collapses and experiences weak accretion shocks. 
However, it is difficult to assess the importance of gradual 
vs. abrupt heating based on the median history of a large 
sample of particles. A better approach is to look at the ther- 
mal history of individual ejected gas particles. 

We randomly select ejected gas particles associated with 
a randomly-selected galaxy group in ZCool.SF.SN. AGN. 
For reference, the randomly-selected group has a total mass 
of M 50 o(z = 0) « 9.2 xlO 13 M Q withT 50 o(.z = 0) « 1.3 xlO 7 
K and S^ooiz = 0) ~ 186 keV cm 2 . (Note that we have also 
randomly selected other groups for comparison and find very 
similar results.). We randomly select 10 particles in several 
radial bins (the same bins as adopted in the right panel of 
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Figure 13. The entropy (top), density (middle), and temperature (bottom) histories of randomly-selected ejected gas particles associated 
with a randomly-selected galaxy group in the ZCooLSF_SN_AGN run. Ten ejected gas particles are randomly-selected in 3 radial bins: 
1 < r/r^oo < 2 (left), 2 < r/rsoo < 3 (center), and 5 < r/r 500 < 10 (right). The thick dashed black curve represents the median entropy 
history of gas constituting the present-day IGrM (i.e., r < rsoo) and the coloured curves represent the 10 randomly-selected ejected 
gas particles. Broadly speaking, based on the entropy histories there are two classes of ejected gas particles: those which were heated 
gradually and those which were heated abruptly. The latter probably correspond to particles that were directly heated by BHs, while 
the former are likely those particles which were heated from collisions with the outflowing directly-heated gas. 



Fig. 11). In Fig. 13 we show the entropy, temperature, and 
density histories of these particles in three of the radial bins 
(see caption). 

Examining first the entropy histories (top panels) we see 
that, broadly speaking, the particles may be broken up into 
two classes: one where the entropy rises gradually, similar to 



(but in excess of) the gas destined to form the IGrM (i.e., 
r < ?"5oo; thick dashed black curve), and another where the 
entropy first drops before being abruptly raised to nearly its 
final value. 

The physical interpretation of the latter population is 
simple: it represents gas that significantly cooled (radia- 
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tively) and collapsed to the centre of its halo before being 
directly heated by a supermassive BH. We note that the 
entropy jump experienced by this gas (often ~ 4 orders of 
magnitude) is about what we would expect given that this 
gas cooled down to ~ 10 4 K (see bottom panels) before 
being heated^ by AT hcat = 10 s K (see Section 2.1). Some 
additional entropy generation occurs after the direct heating 
episode, which is likely due to weak shocks experienced by 
the gas as it flows outward. 

In order to determine the nature of the heating of the 
class of particles whose entropy increases much more grad- 
ually with time (but still in excess of the gas that ends up 
in the group), one needs to also examine the density and 
temperature histories. For this class of particles, most of 
the entropy production occurs between 1 < z < 4. Over this 
time period the density of these particles is typically slightly 
decreasing with time towards the present day. This indicates 
that the gas is outflowing and, as it does so, is experiencing 
numerous weak shocks. In fact, it is likely the gas is outflow- 
ing because it has experienced weak shocks, from collisions 
with rapidly outflowing directly-heated gas. 

That the outflowing, directly-heated gas shock heats 
other gas leading to further ejection can also be arrived 
at through other considerations. In particular, the mass of 
directly- heated gas (niheat) is given b}0 

m hc at = e r e f m B HC 2 /[(3/2)fe B AT heat ] jurorr , (2) 

where e r e/ = 0.015 is the feedback efficiency (see Booth 
& Schaye 2009) and AT heat = 10 8 K sa 8.6/k B keV is the 
energy injected per particle (see Section 2.1). 

Over the relatively small range in halo masses consid- 
ered here, we find that the total mass of BHs within rsoo is 
to-bh(< J'soo) ~ 2 x 10 -5 M5oo- Filling in the other numbers 
yields the approximate scaling between the heated gas mass 
and the group mass: 

m hcat « 1.3 x 10~ 2 Af 50 o (3) 

One can compare this estimate of the directly-heated 
mass of gas with the typical mass of ejected gas (see Fig. 8), 
which approximately scales with M500 as 

m eject ed » A(-0.5Iog 10 [JW B oo/3 x 1O 13 M ] + 0.46) . (4) 

Therefore, at our median mass of M 500 ~ 3 x 1O 13 M 
the ratio of ejected to directly-heated gas mass is « 6.2. This 
indicates that the directly-heated gas is able to heat neigh- 
bouring low-density gas, leading to additional mass outflow. 
But note this process becomes less effective with increasing 

8 The temperature histories typically only show a jump of ~ 
3 orders of magnitude, rather than 4, which would appear to 
be inconsistent with the entropy jump. However, we note that 
the time scale for a single hot particle to cool adiabatically by 
free expansion is ~ h/c s , where h is the SPH kernel length [~ 
(m/p) 1 / 3 ] (where m and p are the gas particle mass and density) 
and Cs the sound speed. For groups with c s of a few hundred 
kilometers per second, this timescale is short compared to the 
interval between snapshots. Thus, we expect directly-heated gas 
particles to cool due to adiabatic expansion. However, this kind 
of expansion does not affect the entropy. 

9 This is actually an upper limit to the mass of directly-heated 
gas, since it neglects BH growth due to mergers and assumes that 
a gas particle can be heated only once. 



group mass. For example, for the randomly-selected group 
studied in Fig. 13, which has a mass of M500 ~ 9.2xl0 13 M o , 
the expected ratio of ejected to directly-heated gas mass is 
« 2.9. This ratio is approximately consistent with the frac- 
tion of particles that experience large entropy jumps in Fig. 
13, confirming our hypothesis that these particles have been 
directly heated by the AGN. 

4.4 Sensitivity to parameters of the AGN 
feedback model 

One of the main motivations for the OWLS project is to 
quantify the effects of varying subgrid physics in the simu- 
lations. This is done by varying the parameters of a given 
subgrid physics module, or by switching on/off the various 
modules one at a time. In a future paper we will present a 
detailed comparison of groups from all of the ~ 50 OWLS 
runs, which, among other things, vary the parameters of the 
prescriptions for stellar populations, chemodynamics, radia- 
tive cooling, and SN and AGN feedback. For the purposes 
of the present study, however, it is worth briefly discussing 
the sensitivity of our results/conclusions to parameter vari- 
ations for our implementation of AGN feedback (see Booth 
& Schaye 2009 for a detailed description of the model). 

As already noted in Section 2.1, a certain fraction of the 
rest mass energy of accreted gas, e, is used to heat a number 
of randomly-selected neighbouring gas particles. Through 
systematic variation, Booth & Schaye (2009) found that a 
value of e ~ 0.015 reproduces the local relations between 
BH mass and stellar mass and velocity dispersion, as well 
as the z = cosmic BH density. We therefore decided to 
leave e fixed at this value. BHs in the simulation store up 
this accreted 'feedback energy' until it is sufficient to raise 
the temperature of riheat gas particles by ATh ea t- We have 
explored the effects of varying these two parameters on the 
properties of the IGrM. In particular, we have tried varying 
ri hcat at fixed AT hcat and Ar hcat at fixed n hcat , as well as 
fixing the product Hheat ATheat (which implies a fixed level 
of energy injection per feedback episode) while varying both 
parameters. 

We find that the results are quite insensitive to varia- 
tions in Wheat at fixed ATheat, but are somewhat sensitive 
to variations in ATheat at fixed ni, eat . This is not unex- 
pected, since if ATheat is set to values that are comparable 
to, or smaller than, the temperature of the local ambient hot 
medium, the heating will be inefficient. This is because the 
gas is heated at fixed density, implying AT/T = AS/S, so 
that a small fractional change in temperature corresponds 
to a small fractional change in entropy. 

In terms of the median entropy of the IGrM, we find 
that the result is insensitive to the exact value of ATh ea t 
so long as it is sufficiently high to expel gas from the low- 
mass progenitors of galaxy groups. As we have shown above, 
it is the ejection of gas from these high-z systems that 
yields the required truncation of the IGrM entropy distri- 
bution. To give one specific example, we found very similar 
median IGrM entropies (compared to ZCooLSF_SN_AGN, 
which uses ATii Ca t = 10 s K and Wheat = 1) f° r a run with 
ATheat = 10 7 K and Wheat = 10. Note, however, that in 
this case ATh ea t is comparable to the virial temperature of 
the z — groups we are studying. The argument presented 
above therefore implies that heating will be ineffective in 
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these haloes today, and indeed we find that the entropy at 
the centers of the groups is lower than observed (even though 
the median entropies of the groups remain similar to what 
is observed). In other words, decreasing Al^eat to this level 
results in a pile-up of low entropy gas, since the entropy 
added to the gas by AGN is not sufficient for the gas to rise 
buoyantly out of the group core. 

The upshot of the above is that as long as the AGN 
is 'able to do something' (i.e., the injected energy is not 
immediately radiated away), it tends to yield the same IGrM 
properties, independent of the exact choice for the values of 
the parameters of the model. 



5 DISCUSSION AND CONCLUSIONS 

In Paper I we demonstrated that the ZCooLSF_SN_AGN run 
that is part of the large OWLS suite of simulations (Schaye 
et al. 2010) simultaneously reproduces an impressively wide 
range of X-ray and optical properties of groups, including 
the entropy, temperature, and metallicity distributions and 
their dependence on halo mass, as well as the stellar mass of 
the group and the central brightest galaxy (CBG) and the 
stellar age of the CBG (see also Sijacki et al. 2007; Puch- 
wein et al. 2008; Fabjan et al. 2010, who also find improved 
agreement with the data when AGN feedback is included 
in their simulations). In the present study we have used a 
number of additional OWLS runs (including the effectively 
non-radiative and no feedback runs NoCool and ZCooLSF, 
respectively) to deduce exactly how the IGrM entropy is 
raised in the ZCooLSF.SN.AGN run. 

Somewhat surprisingly, we find that none of the previ- 
ously proposed mechanisms for raising the entropy of the 
IGrM are a good description of what happens in the sim- 
ulation. The excess entropy of the gas within groups today 
is not the result of heating of that gas. On the contrary, we 
find that the gas within groups today is that which pref- 
erentially avoided strong heating by BH feedback. Instead, 
we find that the excess entropy is the result of selectively 
removing the lowest entropy from the system by star for- 
mation and, much more importantly, by the ejection of gas 
from low-mass progenitor haloes, particularly over the red- 
shift range z ~ 2 — 4 (see Fig. 11). 

The idea of selective removal of low entropy material as 
a mechanism for raising the entropy of the IGrM goes back 
to the studies of VB. Originally, these authors envisaged that 
the lowest entropy material would end up forming stars [al- 
though Voit & Bryan (2001) and Voit et al. (2002) did note 
the possibility that extreme heating could also remove the 
lowest entropy gas]. However, it is now clear that in order 
to reproduce the observed excess entropy via cooling/star 
formation alone, an unacceptably large fraction of the gas 
would need to be converted to stars. Bower et al. (2008; here- 
after B08) proposed a simple modification of the VB trun- 
cation mechanism, in which BH feedback specifically targets 
the lowest entropy gas. In this way, the entropy of the IGrM 
is still raised by the removal of the lowest entropy gas while 
allowing the possibility of retaining reasonable stellar mass 
fractions. Indeed, B08 implemented this simple modification 
of the VB mechanism into the Durham GALFORM semi- 
analytic model of galaxy formation and showed that it can 



yield simultaneous matches to the mean thermal properties 
of the IGrM and to the stellar content of group galaxies. 

While our detailed cosmological hydrodynamic simula- 
tions also show that the primary mechanism for generating 
excess entropy in groups is through ejection of low entropy 
gas, the way in which this happens differs from that pro- 
posed by B08. In particular, B08 (see also Bower et al. 2006) 
implemented a prescription for 'radio mode' AGN feedback 
which becomes effective in high mass systems where quasi- 
hydrostatic gaseous atmospheres can develop. In their model 
ejection only occurs if the energy imparted to the gas exceeds 
the gravitational binding energy of the present-day group. 
Consequently, the energy required to eject the gas is sub- 
stantial (see B08 for further discussion). By contrast, in our 
simulation the ejection occurs at high redshift when the low 
entropy gas is located in the shallower potential wells of the 
low-mass progenitors of groups. Consequently, the energy 
required to eject the gas is lower than that required by the 
model of B08. The ejection in our simulations occurs dur- 
ing the phase when the black hole accretion rate is close to 
Eddington limited. In other words, we find that it is 'quasar 
mode' AGN feedback (as opposed to 'radio mode') that is 
responsible for raising the entropy of the IGrM over and 
above that due to gravitational shock heating. 

In the present study we have of course not explored the 
full range of conceivable AGN feedback models. The simu- 
lations are therefore not guaranteed to capture the physics 
of AGN feedback correctly. However, we are encouraged by 
the fact that the model simultaneously reproduces a wide 
range range of global and structural X-ray/optical proper- 
ties of galaxy groups (see Paper I) with virtually no tuning 
of the model parameters (in terms of AGN feedback, only 
the efficiency of feedback was tuned so that the simulations 
reproduce the normalization of the local black hole scaling 
relations; see Booth & Schaye 2009). A direct test of the 
model will be to search for outflows at high redshift to con- 
firm their existence and validate the model. Indeed, there 
are already signs that such 'quasar mode' outflows are op- 
erating both at low and high redshift (see, e.g., Moe et al. 
2009; Dunn et al. 2010; Tombesi et al. 2010a,b; Feruglio et 
al. 2010; Raimundo et al. 2010). Determining how common 
these outflows are and what their typical properties are (i.e., 
whether they are ejecting sufficient mass on average to ex- 
plain the overall entropy and baryon mass fractions of local 
groups) will be key. 

How realistic is the current implementation of AGN 
feedback given our existing knowledge of real quasar-driven 
outflows? Recent observations suggest outflow velocities up 
to a significant fraction of the speed of light on very small 
scales (~ pc; e.g., Tombesi et al. 2010a). Of course there is 
no hope of being able to resolve such tiny scales in cosmo- 
logical simulations at present. However, the hot, outflowing 
gas will expand and cool adiabatically and as long as: i) the 
simulations capture the outflow when the velocities are still 
supersonic with respect to the the IGrM; and, ii) the cooling 
time of the shocked gas is large; then the simulations should 
retrieve approximately the same result for a fixed amount of 
energy. In this respect, an important parameter of the model 
is AT heat . If AT hoat < T vir , then heating will be inefflcien10 



As we discussed in Section 4.4, in terms of reproducing the 
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(i.e., it violates the second condition), which changes the 
physics and the results will be sensitive to the exact value 
AThcat (see Section 4.4 for further discussion). This is the 
reason why we have adopted a 'high' value for AT hcat - it 
guarantees efficient feedback on the scale of groups (but note 
that this choice in no way guarantees that the simulations 
will reproduce the observations) . 

Another possible way to test the model may be through 
observations of cosmic rays (CRs), since the AGN feedback 
results in some very strong shocks, with Mach numbers ap- 
proaching ~ 100 (see Fig. 13). This may result in the pro- 
duction of CRs, but whether or not this will be easily distin- 
guishable from CRs produced in shocks driven by SNe and 
mergers (e.g., Pfrommer et al. 2007) will require further de- 
tailed modelling. 

It is worth briefly commenting on the suitability of the 
SPH formalism to the problem at hand. Since SPH by its 
nature will smooth out shocks, can we be confident that 
the rate of shock heating in the simulations is correct? We 
have verified that a version of the AGN model run with 
8 times lower mass resolution (2 times lower spatial reso- 
lution) yields remarkably similar entropy distributions for 
the IGrM. Also, as long as the radiative cooling time in the 
post-shock gas is long, the correct result should be obtained, 
since the energy is conserved and the gas will expand in both 
cases (i.e., in a real shock and in a SPH shock). The imple- 
mentation of AGN feedback in our simulations guarantees 
that the post-shock gas has a long cooling time. 

It is interesting to ask the question whether only AGN 
feedback is capable of 'raising' the entropy of the IGrM 
to the observed level. For example, perhaps an alternative 
implementation of SN feedback in the simulations could 
achieve a similar result. We find that it is difficult to rule out 
this possibility based on energy arguments alone. We have 
calculated the total energy injected by BHs within rsoo in 
our simulation as [e r e//(l — e r )]?riBH(< 7*500 )c 2 , where e r 
and e/ are the radiative efficiency of the BH accretion disk 
and the fraction of emitted energy assumed to couple to 
the gas and triBE(< fsoo) is the total mass of BHs within 
rsoo. This yields a typical total injection energy of a few 
10 61 ergs, or approximately 1-2 keV per particle, assuming 
-Mbaryon^soo) = faMzoo- The maximum amount of energy 
available from SNe is estimated as follows. We sum the ini- 
tial (i.e., ignoring stellar mass loss) stellar mass of star par- 
ticles within r 5 oo in ZCooLSF_SN_AGN (which reproduces 
the z = stellar mass fractions of observed groups, see Pa- 
per I). Assuming a Chabrier IMF, we calculate the number 
of stars in the mass range 8 — 100M Q , all of which are as- 
sumed to yield core-collapse SNe. Finally, the energy per SN 
is assumed to be 10 51 ergs. Typically, we find that the en- 
ergy available from SNe is 0.5 — 1.0 times that injected by 
BHs in our simulation (note, however, that the maximum 
amount of energy that can be injected by BHs is an order 
of magnitude larger, but that the adopted BH efficiency was 
tuned to match the normalisation of the z — BH scaling 
relations - see Booth & Schaye 2009). Similar results are 



median entropies of the IGrM, what is impor- tant, is the depth 
of the potential well (characterised by Tvir) of the progenitors of 
groups, rather than that of present-day groups. 



obtained if we use r2oo as our limiting radius, rather than 

T500. 

Energetically, therefore, the contribution from SNe is 
significant. However, whether or not SNe can significantly 
raise the entropy of the IGrM depends on the importance of 
radiative losses. Since the energy injected by SNe will likely 
be much more spatially-distributed than that by BHs, one 
expects radiative losses to be more important for the former. 
Indeed, this is the reason why ZC00LSFJSN yields much 
higher stellar mass fractions and lower IGrM entropies than 
ZC00LSFSN-AGN in the present study. We also note that 
none of the many variations of the SN-driven feedback model 
parameters explored as part of the OWLS project yield group 
properties that are similar to ZCooLSF_SN_AGN (this will 
be presented in detail in a future paper). A wide range of 
implementations of SN-driven feedback have been explored 
in previous studies and similar conclusions have been drawn 
(e.g., Thacker & Couchman 2000; Borgani et al. 2002, 2004; 
Kay et al. 2004). 

Two notable exceptions are Kay (2004) and Dave et al. 
(2008). Dave et al. (2008) implemented a stellar feedback 
model in which the initial wind mass loading and velocity 
scale with halo properties as expected on larger scales for 
outflows driven by radiation pressure from massive stars on 
dust grains. They found groups with stellar mass fractions 
that are similar to the observed mass fractions. It is worth 
noting, however, that the normalisation of their feedback 
prescription was treated as a free parameter and that the 
amounts of energy and momentum injected in their out- 
flows both far exceed those available from star formation. 
Kay (2004) showed that if the energy injected by SNe is 
highly-targeted (i.e., deposited in a very small amount of 
gas) and raises the entropy of the gas to a very high level 
(he adopted 1000 keV cm 2 ), it can yield group properties in 
reasonable agreement with observations. This result is inter- 
esting because in this case the SNe 'behave' much like AGN. 
The downside of this approach is that one must fine-tune 
the parameters of the SN feedback model to achieve a result 
that comes about naturally from a model that includes BH 
growth and AGN feedback. 

In the present study we have concerned ourselves with 
the establishment of only the global (median) properties of 
the IGrM and found that quasar mode AGN feedback is the 
dominant mechanism for raising the overall entropy of the 
gas. However, this statement does not imply that SN feed- 
back or radio mode AGN feedback are unimportant. As we 
have already discussed in Section 3.3, SN feedback prevents 
star formation in high-z galaxies. This gas ends up as part of 
the IGrM but it has not been heated to such an extent that 
it modifies the median entropy of the hot gas. Low accretion 
rate (i.e., radio mode) AGN feedback is likely important 
at late times to prevent the development of cooling flows. 
But it should be kept in mind that by mass the gas within 
the cooling radius (where the radio mode heating tends to 
be targeted) represents only about 10-20% of the total. It 
is therefore difficult for radio mode feedback to affect the 
global properties of the gas, which is necessary in order to 
match observations of groups [see also the recent study of 
Fontanot et al. (2010) who find that semi-analytic models 
that invoke radio mode only, and whose parameters have 
been adjusted to reproduce the galaxy population, tend to 
greatly overproduce the observed radio emission] . 
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Finally, we have focused on the global thermal proper- 
ties of the hot gas in groups and low-mass clusters. However, 
more massive galaxy clusters also show evidence for excess 
entropy, at least in the central regions (see, e.g., Pratt et 
al. 2010). At present, it is widely believed that this excess 
entropy is due to AGN radio mode feedback [in 'cool core' 
(CC) clusters] or cluster-cluster mergers [in 'non-cool core' 
(NCC) clusters]. However, there are problems with both of 
these explanations. In the case of radio mode feedback in 
CC clusters, while there may be sufficient energy to offset 
radiative cooling losses (e.g., Dunn & Fabian 2008), there 
does not appear to be sufficient energy to explain how these 
systems obtained their high entropies in the first place if 
the heating occurred via radio mode after cluster formation 
(McCarthy et al. 2008; NCC clusters cannot be explained 
in this way either). As for mergers being an explanation for 
NCC clusters, Poole et al. (2006, 2008) used high-resolution 
SPH simulations of cluster mergers to show that the dense, 
low entropy gas in the cores of clusters is generally robust to 
mergers (although see Mitchell et al. 2009 and ZuHone 2010 
who argued that SPH simulations may underestimate the 
degree of mixing during mergers). In any case, the physical 
origin of the high entropy in clusters (and its scatter) is still 
very much an open question. 

AGN feedback at high redshift (in quasar mode) may 
provide an alternative explanation for the excess entropy of 
the ICM in both CC and NCC clusters. Naively one might 
rule the quasar mode hypothesis out by noting that BHs are 
not sufficiently energetic to eject gas from the most mas- 
sive systems (which our simulations confirm; see Fig. 3). 
However, we note that, unlike groups, the excess entropy in 
clusters is confined to relatively small radii (the gas at in- 
termediate/large radii is described well by the self-similar 
model). Thus, to explain the excess entropy in clusters the 
gas may not have to be ejected to beyond the virial radius. 
Heating it at high redshift so that it ends up at intermedi- 
ate/large radii may be sufficient. We intend to present the 
implications of our model for massive galaxy clusters in a 
future study. 
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APPENDIX: EFFECTS OF UV/X-RAY 
PHOTOHEATING ON THE THERMAL 
HISTORY OF GALAXY GROUPS 

All of the OWLS runs include photoheating from the Haardt 
& Madau (2001) model for the UV/X-ray background from 
galaxies and quasars. In order to assess the effects of ra- 
diative cooling, star formation, and feedback on the IGrM 
we therefore also included photoheating in our baseline No- 
Cool model. Here we compare the NoCool model to a non- 
radiative simulation which omits UV/X-ray photoheating. 

Fig. 14 compares the thermal histories of the NoCool 
and non-radiative runs. While it is clearly important at early 
times (z > 2), photoheating has no important consequences 
for the thermodynamic properties of the present-day IGrM. 

This paper has been typeset from a TjnjX/ ETjjX file prepared 
by the author. 
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Figure 14. Comparison of the thermal histories of the NoCool 
run which, like all other runs investigated in this study includes 
UV/X-ray photoheating, and a non-radiative run, which neglects 
photoheating. The curves have the same meaning as in Fig. 5. 
UV/X-ray photoheating is unimportant at late times, after group 
formation. 



